Lower Bounds on the Exchange-Correlation Energy in Reduced Dimensions 
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Bounds on the exchange-correlation energy of many-electron systems are derived and tested. By 
using universal scaling properties of the electron-electron interaction, we obtain the exponent of the 
bounds in three, two, one, and quasi-one dimensions. From the properties of the electron gas in the 
dilute regime, the tightest estimate to date is given for the numerical prefactor of the bound, which 
is crucial in practical applications. Numerical tests on various low-dimensional systems are in line 
with the bounds obtained, and give evidence of an interesting dimensional crossover between two 
and one dimensions. 
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In 1979 Lieb [l[ planted a landmark in quantum many- 
body physics by proving the existence of a lower bound 
on the indirect part of the Coulomb interaction. The 
existence of such a bound is of immediate relevance to 
such fundamental questions as the stability of matter Q • 
For the purpose of quantitative calculations, on the other 
hand, existence of a bound is not enough - one would 
wish it to be as tight as possible. A tighter version of 
Lieb's bound was later derived by Lieb and Oxford 
and it is this tighter form, known as Lieb-Oxford (LO) 
bound, which is used as key constraint in the construc- 
tion of many modern density functionals 0, Q ^ which in 
turn are used in calculations of the electronic structure 
of atoms, molecules, nanoscale systems, and solids. 

In connection with recent advances in low- dimensional 
physics it is important to ask whether LO-like bounds 
exist and can be formulated also in reduced dimensions, 
in particular since the study of low-dimensional systems 
today forms a significant part of condensed-matter and 
materials physics. 

The LO bound [3], in its original form, applies 
to all three-dimensional (3D) nonrelativistic, Coulomb- 
interacting systems. The bound can be expressed in 
terms of the indirect part of the interaction energy 

(EI!, 

W X M = h\V e ,\f) - U[n] > - C 3 / dW/ 3 (r), (1) 



where the electron-electron (e-e) interaction operator is 
Coulombic, i.e., V ee — J2i>j l r » ~ r j'l -1, expectation 
value is calculated over any normalized many-body wave- 
function *ff(ri, ...,rjy). n(r) is the corresponding density, 
and U[n] is the classical Hartree energy. For the prefac- 
tor C3 , where the subscript denotes the number of dimen- 
sions D, Lieb originally found C3 = 8.52, which was sub- 
sequently refined by Lieb and Oxford to C3 = 1.68, and 
later, numerically, by Chan and Handy to C3 = 1.64 



0] • Recent numerical studies [H, 0] , as well as modeling 
of the prefactor based on its known properties [l(J , have 
given evidence that the bound can be further tightened. 

In two-dimensions (2D), Lieb, Solojev and Yngva- 
son [ll| (LSY) showed that 



W XC [V] >-C 2 d 2 rn 3/2 (r), 



(2) 



where C2 SY 



192V27T w 481 > C£ u . For a 
dimensional system, the bound may be written as 



D- 



W XC [V] >-C D 



d D r n c 



(r), 



(3) 



but we note that the existence of a bound of this form has 
been rigorously proven for only 3D and 2D, and that the 
tightest possible form (i.e., the smallest possible value of 
C^) is unknown in all dimensions. 

In this paper we (i) show that the exponents of n in 
Eqs. ([1]) and J2} are consequences of universal scaling 
properties of the e-e interaction; (ii) use this result to de- 
duce the exponent a of a possible one-dimensional (ID) 
bound; (iii) provide an estimate of the prefactor Cd that 
corresponds to a dramatic tightening of C^ SY , smaller 
but still significant tightening of C3 , and the first pro- 
posal for C\\ (iv) observe unexpected parameter inde- 
pendence and generality of the bound with respect to 
the model chosen for interactions in ID; and (v) test the 
ID and 2D bounds against analytical and near-exact nu- 
merical data for various low-dimensional systems. 

The ID case, in fact, is subtle because the Coulomb 
interaction is ill-defined. Hence, we consider a contact 
interaction, V ee = vJ2i>j ^( x i ~ x j) with rj > 0. The 
discussion below on the ID case refers to this type of 
interaction. However, we consider also a soft-Coulomb 
interaction, V ee — J2i>j [( x i ~ x j) 2 + ° 2 ] > which corre- 
sponds to a quasi-lD (qlD) situation. 

Under homogeneous scaling of the coordinates, r — * 71- 
(0 < 7 < 00) 6], the (DiV)-dimensional many-body 



2 



wavefunction scales as ^'(ri . . . r^r) — ► v I' 7 (ri . . . tn) = 
ry DN / 2 \$j(jY 1 . 'yrjv), preserving normalization. This 
yields the number-conserving scaled density n(r) — > 
7i 7 (r) = ^ D n("fr). On the other hand, W^d^] — ► 
W xc [* 7 ] = 7W XC [*], since both the Coulomb (D = 2,3) 
and contact (D = 1) interaction, and their Hartree ap- 
proximations scale linearly. Thus, Eq. ([3]) becomes 



(4) 



and consistency between Eqs. ([3]) and (|4]) immediately 
yields 7 = ^ D ( a ~ x \ giving a = 1/D + 1. For D = 3 
and D = 2 this yields a = 4/3 and a = 3/2, respec- 
tively, in agreement with the LO and LSY bounds. We 
thus find that if a bound of this form exists, its expo- 
nent is, in all dimensions, uniquely determined by coor- 
dinate scaling, without requiring the complicated analy- 
sis performed in Refs. 0, i, and [H|. For a LO-like 
bound in ID, the same scaling argument suggests the 
form Wxcf'J'] > — G\ Jd 1 rn 2 (r), although the existence 
of such a bound in ID is at present only a conjecture. 

The exponent a = X/D + 1 in Eq. ((4| is the same 
as in the expression of the exchange energy, E x [n], of 
the homogeneous D-dimensional electron gas, which is 
applied in the local-density approximation (LDA) [12| to 
the inhomogeneous case. Thus, we can express the right- 
hand side of all LO-like bounds in terms of 



E. 



LDA r 



A, 



d D rn a (r), 



(5) 



where A* = Z^n' 1 / 3 /4, A 2 = 2 5 / 2 tt- 1 / 2 /3, and A x = 
r]/4 12]. For ground-state densities, the left-hand- 
side of Eq. (J3j> can be written in terms of the full 
exchange-correlation energy, E xc [n] = W xc [n] +T c [n] > 
W xc [n], where the inequality follows from the positive- 
ness of the kinetic-energy part T c of the correlation 
energy, and the density functional W^n] is obtained 
by evaluating W^c^] with the wavefunction minimizing 



T + V e( 



under the constraint of reproducing the 

ground-state density n(r). We can now cast the bound 
in the form 



E, 



s [n] > \ D E™ A [n] 



(6) 



where Xd = Cd/Ad- Thus, the tightest possible bound 
can be obtained by looking for the maximum value of the 
density functional 



Ad M 



E r 



E x [n] 



E c [n] 



£LDA[ n ] E% DA [n] E^ A [n] 



(7) 



over all possible D-dimcnsional many-body systems. 
Clearly, this maximization cannot be performed in prac- 
tice. However, we can make an educated guess as to what 
the resulting maximum will be. 



First, we note that for 3D, it was shown rigorously that 
the constant A3 in Eq. ([6]) can be replaced by a monotonic 
function depending on particle number, As(iV), which as- 
signs to all systems with particle number N a common 
value A3(iV) < A3 (TV — > 00) = A3, such that a LO-like 
bound with A (TV) in place of A holds for all systems with 
this N I, [lOl. The values commonly quoted for A3 , A3 , 
and A^ 11 are actually estimates of \s(N — ► 00). To obtain 
the tightest possible universal (iV-independent) bound 
we thus take N — > 00. 

Second, we may expect that E x is always relatively 
close to -E^ DA , as the particle density n is varied over 
different physical systems. Hence, the functional Ar)[n] 
is expected to be largest for systems where the right- 
most ratio in Eq. ([7]) is largest, i.e., where correlation is 
largest relative to exchange. This situation is typical of 
the extreme low-density limit n — > 0. 

Taken together, the rigorous property of a maximum 
at N — > 00, and the nonrigorous but reasonable re- 
quirement that n — > 0, suggest that the largest possi- 
ble value of A3I71] is obtained for the r 3D — > 00 limit of 
the 3D electron gas (3DEG) with the density parame- 
ter r 3D = 3 1/3 (4t™)- 1 ' /3 This expectation leads to 
A3 = A 3D EG[ri D -h. 00] = 1 + e c (r 3D -> oc)/e x (r 3D 
00) = 1.9555 [8(, where e x (e c ) is the exchange (correla- 
tion) energy per electron. We note that the original LO 
bound with C%° = 1.68 corresponds to A^ = 2.27, while 
the present estimate A3 = 1.9555 is tighter, and consis- 
tent with the empiricalprefactor obtained by evaluating 
As[n] for real systems [8J. 

We now assume that the above argument about the 
maximum of A3 carries over to reduced dimensions. For 
the 2D electron gas (2DEG) with r 2D = Xj^/wh we have 
e x (rf>) = -4V2/(37rr 2 °L From the Madelung energy 
of the Wigner crystal [l4j we extract that the leading 
contribution to the correlation in the low-density limit 
is e c (r 2D -f 00) = -0.509/r 2D + 0.815/ (r s 2D ) 3/2 . Thus 
we find A2 = A2degHP ~~ > °°] — 1-84, which is a dra- 
matic improvement on the rigorous mathematical result 
of Ref. [lH that A 2 < A^ SY w 452. 

Analogously, for the ID electron gas (1DEG) with a 
contact interaction and r* D = l/2n we find e x (rl D ) = 
— n/(8rl B ). The leading contribution to the corre- 



,1D 



energy in the low-density limit is e c (?-g 
-77/(87-1°) + vr 2 /32(r s 1D ) 2 [15], yielding Ai 



lation 

00) = 

AiDEG[ r s D ~* 00 ] = 2- Note that this is a general result, 
in the sense that it does not depend on the parameter 77 
of the contact interaction. 

Finally, in the qlD system with a soft-Coulomb inter- 
action we first note that the scaling of W xc is nontrivial 
due additional scaling of the parameter a: ^"Jf] — > 
W£ c [%] = jW^}. Second, for the LDA in Eq. © we 
have now A q iD — 1/2, a — 2, and in the limit a — > 
the integrand is multiplied by a function f a [n(x)\ = 
ln[2/7ran(x)] + 3/2 — /j, where fi ~ 0.577 is Euler's 
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TABLE I: Estimated prefactors \d = Cd/Ad for bounds on 
E xc and Wxc, compared to literature values, in different di- 
mensions. The quasi-one-dimensional (ql) case involves spe- 
cial conditions (see text). 
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2 1 
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A D 


1.96 


1.84 2.0 


2.0 


\lit 
A D 


2.27 [3] 


452 [11] 





constant [16| . Note the scaling property f a [n(x)] — > 
f a [n 1 (x)] — 7/ 7a [n(a;)]. Now, we assume that the LO 
bound in qlD also has this form, i.e., the integrand of 
the qlD expression for W xc [Eq. (J3}] is multiplied by the 
same factor f a [n(x)]. Under this assumption, we may 
search for the maximum values for A q iD[^] in Eq. ([7]) in 
a similar fashion as in 3D, 2D, and ID. 

In the low-density limit of the qlD electron gas 
(qlDEG) we have e x {rf D -> oo) = n[ln(a7m/2) - 3/2j 
fi]/2 and e c {rf D -> oo) 



n[ln(an/27r) + 3/2 +/i]/2 0. 
A q iDE G [^ 1D -» oo] = 2. We note, 



Hence, we find A q iD 
again, that the leading contribution to A q iDEG is inde- 
pendent of the softening parameter a of the qlD model. 

Moreover, we note the highly nontrivial fact that the 
leading contribution to Aideg and A q iDEG is the same. 
This encourages us to propose Ad = 2 as the tightest 
general bound for both ID and qlD. It is also impor- 
tant to note that in 2D, ID, and qlD, the correction to 
the leading term is negative, decreasing the value of the 
corresponding Xo for finite (non-vanishing) densities, in 
line with our proposal to extract the maximum of Au[n] 
from the low-density limit of the electron gas. The re- 
sults for the tightest bounds in different dimensions are 
summarized in Tabic [T] 

Next, we test our bounds against analytical and near- 
exact numerical data obtained independently for low- 
dimensional systems, in a similar spirit as was done for 
3D systems in Ref. [8]. In particular, we consider 2D 
parabolic (harmonically confined) and hard-wall square 
quantum dots d/| (QDs), where the density parameter 
can be estimated as r s = iV -1 / 6 ^ 2 / 3 (Ref. 0) and 
r s = (ttA0~ 1/2 £ (Ref. 0), respectively (2D indication 
omitted for clarity) . Here u> is the harmonic confinement 
strength, L is the side length of the square QD, and N 
is the number of electrons. 

Figure QJa) shows \2\rt\ for a 2D Hooke's atom, which 
is equivalent to a parabolic QD with N = 2. Here we 
focus on some of the analytical two-electron solutions in 
the range r s — 0.9 ... 52 derived by Taut The maxi- 
mum value max(A2,Hooke) ^ 1-55 is relatively close to the 
corresponding 3D result max(A3 i H 00 kc) ~ 1-49 (Ref. @]). 
Detailed analysis of the low-density behavior of A2 is 
given below. In the noninteracting (high-density) limit, 
where the correlation energy is zero, we find numerically 
A 2,Ho C Skc = E x/ E x DA ~ L1 ° which is slightly below the 
corresponding 3D value [§] of 1.17. 



(a) n = 2 



r = 52, 11, 
s 2.9, 0.9 



2D Hooke's atom 



(b) N = 6 



r =8.6,6.1,4.7, 
s 3.9,1.9 



r ^0 



parabolic 
quantum dot 



(C) 



■N = 2 

-N = 8 
-N = 16 



r =1.3,0.62,0.44 



square 
quantum dot 



(d) 



2D electron gas 



FIG. 1: (Color online) Values for Aa[n] in different finite and 
infinite 2D systems compared to our limit A2 = 1.84 (thick 
dashed line). 



In larger QD systems (N > 2), most reference data are 
given only in terms of ground-state total energies i?tot> 
whereas the calculation of Arj[n] requires knowledge of 
the exact exchange-correlation energy and the electron 
density n. The exact DFT correlation energy can be com- 
puted as £ c [n cxact ] = £tot[n cxact ] - £^ x [n cxact ] where 
Etot is the exact total energy and EXX refers to exact 
exchange. To estimate Az?[n cxact ], we may then perform 
a self-consistent EXX calculation and calculate 

pEXXr„EXXi 1 77. r„cxacti pEXXr„EXXi 
x _ E x [n j + Etot [n j - E tot [n j 

*D T^T.r.Ar._EXXl ' 



In this work we have performed the EXX calculations in 
the Krieger-Li-Iafrate (KLI) approximation [21| w ithin 
the octopus real-space density- functional code 22J. We 
note that according to our numerical test for the 2D 
Hooke's atom, the estimate in Eq. ((5|) yields generally 
larger values for A d [n] than the definition in Eq. ([7]) . 

Figure QJb) shows results for a parabolic QD with 
N = 6. Here we use as the reference data the varia- 
tional quantum Monte Carlo (QMC) total energies in the 
weak-confinement regime [23j. In Fig. [TJc) we present 
results for square-well QDs with L = tt and varying N. 
Again, we use variational QMC data for the total ener- 



gies [24(. Comparison with fixed r s to parabolic - and 
also to circular-well QDs (not shown) - reveals that defor- 
mation from the circular geometry decreases A2 [n] . Sim- 
ilar decrease in A2 [n] is found if the circular confinement 
is made elliptic in the r s = limit (not shown). 

The r s = limit allows testing also within circular 
confinement by varying the curvature, i.e., the exponent 
in Circular (r) — \r\ a . Interestingly, the largest value for 
A2 is obtained at the smallest a we can numerically con- 
sider, i.e., at a = 0.5, which gives \2[n] = 1.110. Overall, 
the numerical results summarized in Fig. [T] show that in 
both finite and infinite 2D systems, values obtained for 
\2[n] are consistently below our limit A2 = 1.84 (thick 
dashed line). 

Finally we consider the low-density limit of the 2D 
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FIG. 2: (Color online) Values for \2[n] in a two-dimensional 
Hooke's atom as a function of r a = 2~ 1 ^lu~ 2 ^ 3 . The circles 
connected by solid lines correspond to \2[n] in the defini- 
tion of Eq. ([7]) . The dashed line shows the auxiliary quantity 
A*[n] (see text). The dotted line corresponds to the quasi- 
one-dimensional result. 



Hooke's atom in detail. In Fig. [5] we show \2[n] (solid 
line) up to the extreme low-density regime. As expected, 
\2[n] first increases as a function of r s . However, at 
r s ~ 50, we find an local maximum of \2[n] ~ 1.55 
[see also Fig. QJa)] followed by a decrease at higher r s . 
By contrast, if the (2D) LDA exchange in Eq. ((7|) is re- 
placed by exact exchange, \2[n] — > X*[n] := E xc /E x = 
1 + E c /E x , the behavior is monotonic (dashed line) as 
expected. Note, however, that A* [n] is not the quantity 
used in the LO bound, but is used here as an auxiliary 
quantity. Examination of the total electron density in 
the low-density regime suggests that the unexpected be- 
havior of A2 [n] in Fig. [2] is due to the breakdown of the 
2D-LDA. Namely, as the confinement is made weaker, 
the electrons are pushed further apart from each other 
leading to a qlD ring-shaped total density. In fact, the 
low-density regime in Fig. [2] shows reasonable agreement 
between A* and the qlD result (dotted line) deduced from 
Fogler [16] with the parameter a estimated from the low- 
density ring-like model by Taut ■ Hence, it is evident 
that decreasing the density in a 2D Hooke's atom leads 
to a dimensional crossover. 

To summarize, we have shown that the exponents 
in Eqs. {TJ) and ([2]) are consequences of universal scal- 
ing properties of the electron-electron interaction. Wc 
have thus been able to deduce the exponent a of a 
one-dimensional bound. Furthermore, we have provided 
a tightening of the prefactor of the three-dimensional 
bound, a dramatic tightening of the prefactor in two- 
dimensions, and the first proposal for the prefactor in 
one dimension. Unexpected generality of the bound with 
respect to the type of interactions in one- and quasi-one- 
dimensional systems was observed. Our numerical tests 
for low-dimensional model systems are consistent with 



the derivations, all showing \2[n] < A2, and display an 
interesting dimensional crossover in the low-density limit. 
Besides their general relevance in quantum many-body 
physics, these results provide constraints for accurate ap- 
proximations of the exchange-correlation functionals in 
any dimension. 
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